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Abstract 

This article investigates the long-time behaviour of parabolic scalar conservation laws of 
the type dtu + AiVyA{y,u) — AyU — 0, where y € and the flux A is periodic in y. More 
specifically, we consider the case when the initial data is an disturbance of a stationary 
periodic solution. We show, under polynomial growth assumptions on the flux, that the 
difference between u and the stationary solution vanishes for large times in norm. The 
proof uses a self-similar change of variables which is well-suited for the analysis of the long 
time behaviour of parabolic equations. Then, convergence in self-similar variables follows from 
arguments from dynamical systems theory. One crucial point is to obtain compactness in 
on the family of rescaled solutions; this is achieved by deriving uniform bounds in weighted 

spaces. 
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1 Introduction 

The goal of this article is to study the long time limit of solutions of the equation 

dtu + diYyA{y,u) - AyU^O, i > 0, ?/ G M^, (1) 

where the flux A : x E — > is assumed to be T^-periodic with respect to its first variable. 
Here and in the rest of the article, T''^ denotes the A^-dimensional torus, i.e. = (R/Z)^. 

Classical results on scalar conservation laws (see for instance fHl [15]) ensure that the semi- 
group associated with equation ^ is well-defined in L^(R^) + L°^{R^). The case when the initial 
data belongs to U{y) -|-L^(R) (when iV = 1), where U is a viscous shock profile of equation ^ has 
already been dealt with in a previous article, see [7J. In the present paper, we restrict our study 
to the case when the inital data belongs to v{y) + L^(M.^), where w is a given periodic stationary 
solution of ([l]). 

When the fiux A is linear, say 

A{y,u) = a{y)u, 

this study coincides (at least for some particular functions a) with the one led by Adrien Blanchet, 
Jean Dolbeault and Michal Kowalczyk in [4] on the large time behaviour of Brownian ratchets, as 
we will explain in Remark [2l It is proved in [4] that if the fiux A is linear and if 

limsup ^ / \u{t, y) -v{y)\ {y - ct)^ dy < oo (2) 
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for some velocity c £ which will be defined later on (see ([7])), then there exists a constant C 
and a number k G (0, 1/2) such that 



u{t,v) - v{y) 



Mfoiy) _pf y-ct 



dy < Ct-", 



(3) 



where /o is the solution of an eUiptic equation in (see ([6])), -F is a Gaussian profile, and M is 
the mass of the initial disturbance, i.e. 



M : 



Unfortunately, as we explained in [7j, the above result does not imply that the same convergence 
holds in the nonlinear case. Moreover, the proof of [4], which is based on entropy dissipation 
methods together with Log-Sobolev Poincare inequalities, can hardly be transposed as such to a 
nonhnear setting, although attempts in this direction have been made, see for instance [lO]. Hence 
we have chosen here a slightly different approach, which enables us to extend some of the results 
of ^ to a nonlinear context. Additionally, we recover a weaker version of the convergence but 
without the need for assumption In fact, we prove that ([2]) holds for a large class of initial 
data in v + L^{R^). 

The present work is also embedded in the broader study of the long time behaviour of conser- 
vation laws. We refer the interested reader to the review paper by D. Serre [25 (and the references 
therein) for a thorough description of the homogeneous case, in which the author investigates 
the stability of stationary solutions of scalar conservation laws in various models (parabolic and 
hyperbolic settings, relaxation models...) 

Before stating the main results of this paper, let us now recall a few properties of equation ([T]). 
First, according to a result of |8j, periodic stationary solutions of H]) exist, provided the flux A 
satisfies some growth assumptions. In fact, several different growth regimes were studied in [5]; we 
only recall one of them here, which is the most relevant with regards to our purposes. In the rest 
of the article, we assume that A belongs to Wj];^{T^ x E)^, and that 



3po eR,yye T^, divyA{y,po) = 0. 



(4) 



We also assume that there exists n e (0, {N + 2)/N) such that 



VP > 0, 3Cp > 0, V(p, q) e M^ \p\ < P, 

r \dj,A{y,p + q)-dpA{y,p)\<Cp{\q\ + \qn, 

\ \dWyA{y,p + q)-dWyAiy,p)\<Cp{\q\ + \qn. 

These assumptions were introduced in [71 [8]. They ensure that for any (7 e M, there exists a unique 
periodic stationary solution of ^ with mean value q; we refer to [8] for a discussion of the optimality 
of conditions Moreover, if m is a solution of ([!]) with initial data uu—q 

where v € W^-'°°{R^) is any stationary solution of then u e L°°i[0, 00) x R^). This result will 
be used several times in the article, and its proof is recalled in Appendix A. 

We now introduce the profiles which characterize the asymptotic behaviour of the function u; 
first, the function /o occurring in |[3]) is the unique solution in H^{T^) of the equation 

- Ayfo + div.y(ai/o) = 0, (/o) = 1, (6) 

where ai{y) := {dpA){y,v{y)) G i°°(T^)^. Above and in the rest of the article, the notation (•) 
stands for the average on the torus T^, that is 
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The drift velocity c is then defined by 

c = 7V(ai/o). (7) 

The last function which will appear in the asymptotic profile of u is the equivalent, in the non- linear 
case, of the Gaussian profile F occurring in ([3]); it is the unique solution, in a suitable functional 
space, of an elliptic equation of the form 

- V ?7,j5j9jFm - div^(a;fM) + a • V:ri^M = in R^, with / Fm = M e R, 

l<.,,<iV ■ -^K" 

where the coefficients rjij and a are constant, and the matrix {r]ij)i<ij<N is coercive. Unfortu- 
nately, giving the precise definition of rjij and a would take us too far at this stage. We merely 
recall that thanks to a result of J. Aguirre, M. Escobedo, and E. Zuazua (see flj), the above 
equation has a unique solution for all M e M, and we refer to the next section for more details. 

The main result of this paper is the following: 

Theorem 1. Let A e Wi^^{T^ x M)^, and assume that A satisfies g]), ^. 

Let V be a periodic stationary solution of and let Umi G v + i^(M^). Let u he the unique 
solution of ([T]) with initial condition U|(_q = Umi- Set 

M -.^ I {uini - v) dy. 

Then as t —>■ oo, 

y-ct 



u{t,y) - v{y) - o-,\N/2 -My)Fi 



dy^O. 



Remark 1. In fact, the regularity assumptions on the fiux A are not as stringent as stated in the 
Theorem above. In particular, the conditions on the derivatives with respect to the space variable y 
can be considerably reduced. When looking closely at the proof, the correct regularity assumptions 
on A are 

d^Ae L^^{T^ xR)^ Vfce {0,1,--- ,4}, 
dWyA,divyd^A e i,^,(T^ X M). 

Remark 2. Let us now make precise the link between brownian ratchets and equation ^ in the 
Hnear case. In [4], A. Blanchet, J. Dolbeault and M. Kowalczyk study the long time behaviour of 
the solution f = f{t, y) of the equation 

dtf - A,/ + div,(VV(y - wt)/), i > 0, y e M^, (8) 

with V e C2(T^), w e R^. Setting 

u{t,y) = f{t,y + ujt) Vi>0, VyeR^, 

we see that u satisfies 

dtu + divj,(a(y)u) — AyU = 0, 
where the drift coefficient a is given by 

a(y) = -LJ-Wy^j{y). (9) 

Hence the study of |[8]) and that of |(T|) in the linear case are closely related; they are strictly 
equivalent in dimension one, since any function a G C^(T) can be decomposed as 

a = / a + ia — / = / ct + dycj), for some (j) G C^(T). 
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The equivalence does not hold when N > 2, but in fact, all the results of [4] remain true for an 
arbitrary drift a £ C^(T^) (using exactly the same techniques as the ones developed in [4J). This 
will be a consequence of the analysis we will perform in the next sections. The choice for a function 
a with the structure ^ stems from physical considerations (see [5]): equation |[8| describes the 
evolution of the density of particles in a travehng potential, moving with constant speed uj. 

In the course of the proof of Theorem [U we will also prove that condition ([2]) holds for a large 
class of initial data. The precise result is the following: 

Proposition 1.1. Assume that the flux A is linear, and that Uini G u + L^{R^) is such that 



Let u he the unique solution of ^ with initial data Ui^i. Then ^ is satisfied. As a consequence 
(see l^), {HI holds. 

Hence for linear fluxes and for a large range of initial data, a rate of convergence can be 
given. The derivation of convergence rates in the non-linear case goes beyond the scope of this 
article; in fact, the standard methods to derive convergence rates rely on the use of entropy-entropy 
dissipation inequalities (see [10] in the case of the Burgers equation) , which we have chosen not to 
use here. 

Another consequence of Theorem [T] is the stability of stationary shock profiles of equation H]) 
(see [7j) in dimension one: a stationary shock profile is a stationary solution of |(T]) with = 1, 
which is asymptotic as y ^ ±oo to periodic stationary solutions of |(T|). It was proved in [7] that 
the stability of shock profiles is a consequence of the stability of periodic stationary solutions. Thus 
we have the following 

Corollary 1.1. Assume that N = 1, and that the hypotheses of Theorem[I\ are satisfied. Let 
U e L°°(M) be a stationary shock profile of Jl]). Let Umi G U + L^(M.) such that 



The strategy of proof of Theorem [T] is close to the one developed in |12], in which M. Escobedo 
and E. Zuazua study the long time behaviour of a homogeneous version of ([l]); we also refer the 
interested reader to [H], in which M. Escobedo, J.L. Vazquez and E. Zuazua extend the analysis 
performed in [12] to the case when the flux has sub-critical growth. The flrst step of the analysis 
consists in a self-similar change of variables, which helps us to focus on the appropriate length 
scales; this will be done in the next section, in which we also derive the equations on the limit 
profiles /o and Fm- Then, in section [51 we obtain some compactness on the rescaled sequence 
by deriving some uniform bounds in weighted spaces. Eventually, we conclude the proof in 
Section [4] by using semi-group arguments inherited from dynamical systems theory. 

Throughout the article, we will use the following notation: if tp £ L^^{M.^), we set, for all 
P G [l,oo). 



3m> N + 8, 





and let u he the unique solution of ((Tj) with initial data Ui. 



Then 



lim \\u{t) - ?7||ii(R) = 0. 




and = ||u|li2(^) + \\^ufl^^^^y 
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Sobolev spaces of the type W^'^{^), H''{ip), with s g N arbitrary and p G [1, oo), are defined in a 
similar fashion. When we write ||m||p, or |1u||lp, without specifying a weight function, we always 
refer to the usual norm in E^, with respect to the Lebesgue measure (i.e. ip = 1). 



2 The homogenized system 

The goal of this section is to analyze the expected asymptotic behaviour of the solution u{t) of 
equation (IT|); to that end, we change the space and time variables and introduce a parabolic scaling, 
which is appropriate for the study of the long time behaviour of diffusion equations. Then, using 
a two-scale Ansatz in space and time which was introduced in [4j, we construct an approximate 
solution of the rescaled system. Eventually, we recall and derive several properties of the Hmit 
system. 



2.1 Parabolic scaling 

Consider the solution u £ i,^^([0,oo) x R^) of with u\t=o = "im G v + L'^ n L°°(R^). It is a 
classical feature of scalar conservation laws that the semi-group associated with ^ is contractant 
in Li(R^). Hence, for all t > 0, u{t) ev + L\R^), and 

\\u{t) - u||i < \\u^n^ - v\\i. 



Thus it is natural to compute the equation satisfied by f{t) — u{t) — v E L^(K^) 
stationary solution of |[T]), there holds 



smce w IS a 



dtf + diyyB{yJ)~Ayf^O, t>0,ye 
where the fiux B is defined by 

B{y, /) = A{y, v{y) + f) - A{y, v{y)), y{y, /) e x 

The flux B{y, /) vanishes at / = 0, for all /. Moreover, if the flux A satisfles the assumptions of 
TheoremlH there exists ai € Ci(T^) and Bi e C{T' 



such that 



B{yJ)^a,{y)f + B,{yJ), 



and the flux Bi is such that 



VX>0, 3Cx >0, V/e VyeT'\ Bi{yJ) < Cx\f\. 

At some point in the proof, we will need a more refined approximation of i? in a neighbourhood of 
/ = 0; we thus also introduce aa, as G L°°(T^), B3 e i°°(T^ x M) such that 

B{y, f) = ai(y)/ + a2{y)f + a^{y)f + Mv, /), 

and the flux B^ is such that for all X > 0, there exists a constant Cx > 0, such that for all 

/ e [-X,X], for all y G T^, 

B3{yJ)\<Cx\f\\ 
divyB^{yJ)\<Cx\f\\ 
dfMyJ)\<Cx\f\'. 

The existence of ai (i — 1,2,3) and the bounds on Bi^B^ are ensured by the assumption that 
^ e W'5'°°(T^ X E). Notice in particular that 



ai{y) = a/B(y,/)|/=o = {dpA){y,v{y)), Vy £ T 



N 
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As explained in [4], the interplay between the diffusion and the drift a\ induces a displacement 
of the center of mass. In the linear case, that is, when B\ — 0, the evolution of the center of mass 
can be computed as follows: since the function / satisfies 



9t/ + div„(ai/)- A,/ = 0, 
yf{t,y)dy^N [ a,iy)fit,y) dy. 



there holds 

d 
dt 

Now, for t > 0,?/ e T^, set 

nt,y)^ f{t,y + k). 

feez" 

Since the function ai is periodic, / satisfies 

dtf + divyiaj) - Ayf = 0, t > 0, y G T^, 

and we have, for alH > 0, 



Using Lemma 1.1 of [17] together with a Poincare inequality on the torus T^, it can be easily 
proved that as i — > oo, f{t) converges with exponential speed in L^(T^) towards ^/^ /o, where /o 
is the unique solution of ([6]). Additionally, notice that 

/)= E (/(• + fc))= / f-M. 

Consequently, setting 

c:=iV(ai/o) 

we infer that in the linear case, 



— I {y — ct)f ^ exponentially fast. 



dt 

In fact, it turns out that the nonlinearity has no effect on this displacement, although this is 
not quite clear if we try to include the quadratic term Bi in the above calculation. We will justify 
this result by formal calculations in the next paragraph. Nonetheless, it can be proved in the case 
iV = 1 (see for instance [7]) that when ||/o||i is not too large, 

/ol 



II/WIIl^w yt>o, 



and more generally, the norm of f(t) vanishes for all p € (l,oo]. This somehow explains 
why the quadratic term does not modify the motion of the center of mass for large times: the 
term _Bi(-, /(i, •)) vanishes in L^(R) as i ^ oo. Hence, hereinafter, we choose to make in the 
general case the same change of variables as the one dictated by the linear case. Precisely, let 
U G L,^^([0,oo) X R^) such that 

y-ct 



fit,y) = ^^^±^^U(lo,VTT2t,^^), t>0,yeR-. (10) 

This change of variables is classical in the study of long-time paraboHc dynamics, see for instance 
|12) . In the present case, our change of variables is exactly the same as in [4j; straightforward 
calculations lead to 

drU - div,{xU) + Rdiv,{ia,{z) - c)U) - A,U = -R^+^div^B, (z, , (11) 
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with T > 0, X £ R" , and where 

~ 1 

R = e^ and z — Rx + c — ^ — ■ 

Studying the long time behaviour of / amounts to studying the long time behaviour of U . Now, 
as r oo, the quantity R becomes very large, and thus the variable z is highly oscillating. 
Hence, as emphasized in [4], the asymptotic study of equation ifTTI) somehow falls into the scope 
of homogenization theory; the small parameter measuring the period of the oscillations is then 
£ = R~^ = e""^. However, one substantial difference with classical homogenization problems is that 
the small parameter depends on time, which sometimes makes the proofs much more technical. 
We refer to [4] for more details. 

Let us also mention that the homogenization of equation (fTTj) with a "fixed" small parameter, 
and when the quadratic flux Bi vanishes, has been performed by Thierry Goudon and Frederic 
Poupaud in [14]. As a consequence, the formal asymptotic expansions which will be performed in 
the next section are in fact very close to the ones of [14] . 

2.2 Formal derivation of the limit system 

As usual in homogenization problems (see [3] for instance), the idea is now to assume that the 
solution U of (flTI) admits an asymptotic development in powers of the small parameter measuring 
the period of the oscillations; in the present case, the small parameter is , so that we expect 
the approximation to be vaHd for large times only. Hence, assume that when t 1, 

U{t, x) w Ua (r, x, z) + e"^C/i (r, x, z) + e^^r^^ (-^^ z) H (12) 

where z — e^x + c(e^'^ — l)/2 stands for the fast variable and where for all (r, x) G R+ x M^, the 
function 

Z ^ Ui{T,X,z) 

is T^-periodic. Plugging the Ansatz (fT2|) into equation ifTTjl and identifying the powers of i? = e'^ 
leads to a cascade of equations on the terms Uo, Ui, etc. Notice that according to Lemma A.l in 
the Appendix, / e L°°([0,oo) x M^), and thus U/R^ is bounded in L°°. 

• Terms of order R^: Identifying the highest order terms in equation l|lip when U is given by 
lfT2|) leads to 

c • V.iJo + div,((ai - c)Uo) - A,Uo = -A.C/q + div(aiC/o) = 0, ze . 

We recall the following result, which is a straightforward consequence of the Krein-Rutman Theo- 
rem (see [9]): 

Lemma 2.1. Let a E L°°{T^). Consider the vector space 

E[a] := {w e iJ^(T^), -A^^w + div^(aw) = O} . 

Then dim£'[a] = 1, and there exists a unique function m G E[a] such that (m) — 1. 
Moreover, m G W^^'^iT^) for all p < 00, and 

inf m > 0. 

zeT« 

In the present case, i?[Q;i] — M/o, where /o is defined by |(6]). Hence there exists a function 
F ^ F{t, x) such that 

Uo{t,x,z) = fo{z)F{T,x) V(r,a;,z) G [0,00) X X T^. (13) 
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• Terms of order R^: Concerning the terms of order — e'^ , the case when the space dimension 
is equal to one has to be treated separately. Indeed, 



i?™+Miv,i?i [z,^] = R''''div4a2U') + R'-'^'dW^asU') 



+i?^+Miv.B3(z,-^ 



and using the bounds on B, 



3, 



^) = R''idiv.B,)(^z,^^+RV,U-iduB,)(^z,^ 
= 0{R-'^'^)+0{R-^'^+^)=0{R-^). 



We infer that if U is given by lfT2|) , 



i?^+Miv,^i (^z, = i?2-^div,(a2C/o') (14) 

+ [div,(a2{/o) + 2div4a2C/oC/i)] (15) 

+ i?2-2^div,(a3[/o^) 
+ 0{R-^). 

Consequently, we obtain that when N >2, the term Ui solves the equation 

N Q2TJ 

- A,Ui + div^aiC/i) = -div,((ai - c)Uo) + (16) 

i—l * * 



Since Uo{t,x,z) — fQ{z)F{t,x), we have 

{{ai-c)Uo{t,x,-)) = F{t,x){{aifo) -c)^0 



by definition of c. Hence the right-hand side of l|T6|) has zero mean value, and the compatibility 
condition is satisfied. Thus for all {t,x) £ [0,oo) x R^, ||16p has a unique solution in i/i(T^). 
Moreover, using the linearity of lfT6|) together with the expression (flSl) . we infer that Ui can be 
written as 

Ui{t,x,z)^ fi{z)-V,F{t,x), (17) 

where /i G ^^(T^)^ satisfies 

-AJl,,+dW,{alh,,)^-fo{al.^-c^) + 2^,Jo, V^ G {1, • • • , TV}. (18) 

Notice that according to the regularity assumptions on the fiux A, the function ai belongs to 
^i.oo(T^); thus fo G VF2'P(T^) for all p < oo, and therefore /i G W^'PiT'^) for all p < oo. In 
particular, /i,/o G VFi'°°(T^). 

K = 1, on the other hand, the corrector Ui solves the equation 

- d,,Ui + d.iaiUi) = -d4iai - c)Uo) + - d,{a2U^). (19) 

Notice that the compatibility condition is satisfied, for the same reason as before. Hence in this 

case, 

Uiit,x,z) = fi{z)d,F{t,x) + gi{z)F{t,xf, (20) 

where gi G H^{T) solves 

-Azgi +dz{aigi) = -dz{a2{fof). 



8 



The fact that the compatibility condition is satisfied in all cases justifies the use of the change 
of variables (fTO|) in the nonlinear case. This means that, at least on a formal level, the displacement 
of the center of mass of the function / is unaffected by the presence of the quadratic term Bi . 

• Terms of order As we identify the terms of order one in equation ifTTj) . we obtain 



+div(aiC/2) (21) 

dxidzi 



d U 

-drUo + dW^{xUo) + A^Uo - div^((ai - c)Ui) + 2^ \ + Anl 



i=l 

where the term Anl stems from the expansion of the nonlinear term Bi. According to (fT4|). we 
have 

Anl = d,{a2U^) + 2d,{a2U^U^) + d,{a^Ul) if TV = 1, 

^iVL = div,(Q2C/2) ifiV = 2, 
Anl =0 if TV > 3. 

The evolution equation for the function F follows from the compatibility condition; precisely, we 
obtain 

drF - div,(a;F) - A,F + div, ((ai - c)Ui) - {Anl) = 0. 

We now distinguish between the cases N >2 and = 1. 

[> If A'" > 2 , (Anl) = 0; using |(T7|), we infer that F satisfies 

drF - div^ixF) - ^ =0, t > 0, a; S with Af > 2, (22) 

i<,:j<A' ox.dxj 

where the coefficients {'r]i.j)i<ij<N are given by 

= - (("1,4 - Cj)/lj) . 

The following Lemma entails that equation (|22l) is well-posed (see also Lemma EH]): 

Lemma 2.2. The matrix rj :— {riij)i<ij<N is coercive. 

Lemma [2?2l is proved in [14] in dimension A^, and its proof is recalled in ^4] when A^ = 1. For the 
reader's convenience, we sketch the main steps of the proof here, and we refer to [14J, Proposition 
4.6 for details. 

Proof. Let L be the differential operator 

L4> = ~Az(j) + divz(ai0). 
The idea is to introduce, for all j G {1, • • • , A^}, the function Xj which solves the adjoint problem 

L*X] = -AzXj - oil ■ ^zXj = "ij - (xj) = 0. 

Since the right-hand side satisfies ((aij — Cj)?/;) — for all i/" G kcrL = M/o, the function Xj is 
well-defined. For all ^ € M''^, we have 



= (x-C(-/o(ai-c)-e + 2V./o-e)) 
= -(/oX-^i*(x-e)) + (2x-^V,/o-^) 

= -(i(/ox-Ox-0-2(M-v.(x-0>- 
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Expanding L{fox ■ and using the identity L/q = leads to 

(i(/ox-Ox-e) = (/o|V,(x-Ol')- 

Hence 

E ^^.^-^^^j- = iei' + (/o|V.(x-Ol') + 2(M-V.(x-e)) 

l<ij'<iV 

2\ 



,/oie + v.(x-OI 

We deduce that 

l<i,i<JV 

Now, let ^ £ such that = 0- Since fo{z) > for all z, we infer that 

Taking the average of the above inequality on leads to ^ = 0. Hence the matrix {mj) is 
coercive. □ 

t> If TV = 1 . we have 

(Anl) = a, (aaC/o') = («2/o) 
Moreover, in this case Ui is given by l(20|) : hence 

(div,((ai - c)!7i)) = ((ai - c)/i) + ((ai - c)5i) 

Consequently, the compatibility condition reads 

drF - d^{xF) + ad^F^ - rjd^^F = 0, r > 0, a: G M, (23) 

where the coefficients a, rj are given by 

a (a2/o) + (("i - c)gi) , 
7? := 1 - ((ai - c)/i) . 



Lemma 1231 states that the diffusion coefficient is positive. 

This completes the formal derivation of an approximate solution. In the following paragraphs, 
we recall or prove several results concerning the well-posedness and the long time behaviour of 
equations l(22|) and l(23|) . We will often refer to the equation on F as the "homogenized equation"; 
this term denotes equation (|22l) when N > 2, and (|23l) when = 1. 



2.3 Existence and uniqueness of stationary solutions 

This paragraph is concerned with the existence and uniqueness (in suitable functional spaces) of 
stationary solutions of the homogenized equations (|22|) and (|23| . In the case when N — 1, or when 
iVi,j)i<i,j<N — A/ for some A > 0, such results are stated in [l]. In the general case, we merely use 
a linear change of variables, and the problem is then reduced to the case of an isotropic diffusion. 

Lemma 2.3. Assume that N > 2. For 7 > 0, set tp^ : x E 1-^ exp{j\x\'^). Then there exists 
7 > such that for all M e M, there exists a unique function Fm G H'^(ijj^) satisfying 



l<i,j<N 

Furthermore, the following properties hold: 



ri^jd^djFM - diY^^xFM) = 0, / Fm 



M. (24) 
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N )' 



dF{x) 
dxi 



Thus, always setting x = Py, 

x-\/^F{x) 



N 



Y,{PyUp-'e.)-yyF{y) 



i=i 



v,F(y) 



(25) 



(i) For all M £ R, Fm = Mhi; 

(ii) hi e W^'P n C°°(R^) for all p G [1, oo), and hi € H^{ipj); 
(Hi) hi{x) > for all x e R^. 

Proof. The idea is to perform an affine change of variables in order to transform the diffusion term 
into a laplacian. Precisely, set 

^^_!M±I!2A^ l<^,J<N. 

Then the matrix S — (sij) is symmetric and positive definite (see Lemma [23| : hence there exists 
an orthogonal matrix O € A^n(R) and positive numbers Xj such that 

S' = 0'^Diag(Ai,--- ,Xn)0. 

Let us change the variables by setting 

X = Py, with P := O^Diag(Ay^ • • • , A 

and for any function F € L^(M^), define 

F{y) = F{Py). 
It can be readily checked that for all x e , 

Y] m,] if 5 (x) = y^Vkj f ^ (P^^x), 
j-' " dxidxj ^ oykdyi 

where the coefficients fjkj are given by 

f|k,^J2^p-%AP-\,^^., = {p-'s{p-Y)k,i. 

Using the definitions of the matrices P and S, we infer that 

fi = p-^s{p-Y ^In- 

Thus the diffusion term is transformed into a laplacian with this change of variables. 
Let us now compute the drift term. We have 

= {Py)i, 

and, denoting by (ei, • • • , bat) the canonical basis of M^, 



= {p-'Py)-WyF{y)^yWyF{y). 
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Notice that this property is in fact independent of the definition of the matrix P. Consequently, 
Fm is a solution of l(24|) if and only if Fm satisfies 

-A,Fm - div^(yFM) = 0, / Fm = (det Sy'^^M. 

The only solutions of the above equation in H^{R^) are Gaussian functions. Hence there exists a 
unique solution of l(24l) in H^{R^) for all M, and this solution is given by 

FM{x)^CMexp (^l\P'^x\^ 

where the positive constant C is a normahzation factor. Moreover, 



\P-'x\' 



]^/2 l/2\ ^ 

Diag ( Ai ' ,Xj^' jOx , 



and thus, since \Ox\'^ = \x 



2 



-1/2 / x-l/2 

max A, ) bp < iP^^xP < min A, ) 

l<i<JV / \l<i<N 



All the properties of the lemma follow, with 

7 < 2 max Ai 
\^i<i<Ar J 

□ 

In the case when A'' = 1, the existence of a stationary solution is treated in jl]. Hence we merely 
recall the main results of [T] in that regard. 

Lemma 2.4 (Aguirre, Escobedo, Zuazua). Let M £ R be arbitrary, and let a £ M., rj > 0. Let 
Then there exists a unique function Fm G H^{ipt) which satisfies 

-ridxxFM - dxixFM) + ad^Fl = 0, / Fm = M. 

Moreover, Fm enjoys the following properties: 

(i) Fm e T4^2,P n C°°(M) for all p G [1, 00), and Fm G H'^i^j); 

(ii) If M > 0, then Fm{x) > for all x eR^. 

We deduce from the above Lemma that if 7' < 7, then there exists a constant Cj> such that 

|-FM(a;)| , \dxFM{x)\ < Cy expi-^'x^) ^x e M. 

Indeed, since £ H^^ip^), it can be easily proved that FmV'7' ^ H^{R) for all 7' < 7. Sobolev 
embeddings then imply that FmiPy e W'^'°°(R). 

The existence of stationary solutions of (|22l) and l(23|) is now ensured. We now tackle the study 
of the properties of equations l(22|) and l|23p , focusing in particular on the long-time behaviour and 
on regularity issues. 

We begin with a definition of the weight function K e C°°(]R^), which plays a central role in 
the theory of existence. We use the change of variables l(25l) . which was introduced in the proof 



12 



of Lemma [231 This allows us to transform the matrix {rjij) into the identity matrix. For t > 0, 
y e M^, set F{t, y) = Fit, Py). If is a solution of then F solves 

dtp ~d\Yy{yF)- ^yP = {). 

Consequently, the results of [12] can be directly applied to P, for which existence is proved in the 
functional space L?{Kq), where KQ{y) = exp(?/^/2). Performing the inverse change of variables, it 
is clear that the relevant weight function is given by 

K{x) -.^ Ko{p-^x) = exp 

Notice that by definition of the matrix P, there exist positive constants 7,7' such that 

exp(7'a;2) < K{x) < exp{jx'^) \/x e M^. 
When N = 1, the weight function K is given by 

exp (^). 

We immediately deduce from [l^ the following Proposition: 

Proposition 2.1. Let Fini G L°^{M.'^) n L'^{K). Then the homogenized problem has a unique 
solution 

F e C([0, 00), L^{K)) n C((0, 00), H^{K)) n C^{{0, 00), L^{K)) 

such that i^f=o = Fini- 
Moreover, 

lim \\F{t) - Pm||li(R") = 0, 

where Fm is the unique stationary solution of the homogenized problem with mass M = J^,^ Fini ■ 

Consequently, the homogenized equations (|23|) and l|22p are well posed. We conclude this section 
by stating a result on the construction of an approximate solution: 

Definition 2.1. Let F G C{[0, 00, L'^{K)) nC{{0, 00), H'^{K)). We define the approximate solution 
of (fTTj) associated with F by 

U''pp[F]{t, x; R) = Ua (r, x, z) + R-^Ui (r, x, z) + R~'^U2 (r, x, z) , 

with T >0, X E R^, R > and z := Rx + , and where 

• Uq is defined by (fT3|) ; 

• Ui is defined by ^ if N > 2 and by ^ if N = 1; 

• U2 is defined by 

-A,U2 + diyiaiU2) = (foiz) - 1) [-drF + dW,{xF) + A,F] 

+ (divj;((ai - c)Ui)) ~ divx{{ai - c)Ui) 

i—1 
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Notice that we do not require, in the above definition, that F is a solution of (|22l) or l(23| : hence 
the right-hand side in the equation on U2 is sHghtly modified, so that the compatibiHty condition 
is satisfied and U2 is well-defined. Of course, if F is a solution of (|23|) or i(22|) . the equation on U2 
becomes l(2T|) . 

We then have the following result: 

Lemma 2.5. 1. Let M G M &e arbitrary. Define the function U G L°^{[0,oo) n K^)) by 

U{t,x) -.^U-ppIFm] {T,x;e^). 

Then U is a solution of 
drU - 6iY^{xU) - A^U + Rdiy^{{ai{z) - c)U) = 

where the remainder term JJ'^^'^ is such that there exist C > 0, 7 > such that 

r'"'''"WllLi(e.M^) + r'"'''"WllL=»(M«) < Ce-^ Vr > 0. 

2. Let F^m e i°°(K^) n L'^{K), and let F e C{[0, 00), L'^{K)) be the unique solution of the 
homogenized equation such that F^t^Q = F^i- Let p € C^{M.^) be a mollifying kernel (p > 0, 
J p — I), and let Fg := F *x ps, where ps — 5^^ p{- / 5) , for 5 > 0. 

Let (t„)„>o be a sequence of positive numbers such that lim r„ = +00. For n e N,(5 > 0, 

n — *oo 

define the function uf^ by 

ui{T,x) = U^PPlFs] (T,x;e^"+^) , x e M^, r > 0. 

Then satisfies, with Rn = eJ"'^'^ and Zn = RnX + c ^'"^ , 

druf^ - div2;(xM^) + i?„div^((ai(z„) - c)w* ) - A^u^ = 

= -i?,r+Miv.Si (^z,^^+ri, 

where the remainder term satisfies, for all T > 0, 

\\ri\\L'^([a,T],L\R^y) < wt((5) + Cs^re^''", 

where lot ■ ^+ ^ R+ is a function depending only on T such that linio+ ujt = 0. 

The proof of the above Lemma follows the calculations of the first paragraph; the proof is 
lengthy but straightforward, and is therefore left to the reader. The fact that C/"^™ has exponential 
decay is a consequence of Lemmas 12.31 12.41 

3 Weighted bounds for the rescaled equation 

As explained in the previous section, we choose to work with the rescaled equation ifTTj) rather 
than with the original one Jl]). In fact, it can be easily checked that Theorem [T] is equivalent to 
the following Proposition: 
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Proposition 3.1. Let Uini € L-'^(M^), and let M :— t^mi- 

Let U S C([0, oo), L^(R^)) be the unique solution of l|lip with initial data C/|r=o = Uini- Then 



lim 



C/(t, x)-fo[ e^x + c^-—^ ] FMix) 



dx = 0, 



where the speed c is defined by and Fm € i^(R^) is the unique stationary solution of the 
homogenized equation l|22p . (|23l) with total mass M. 

In turn, since the function /o e L°°(T^) is such that infjN fo > 0, the above statement is 
equivalent to 

lim \\V{t) - Fm||li(r") = 0, 
where the function V = V{t, x) is defined by 

Vir, x) := r > 0, X e R^. (26) 

fo [e^x + c^^) 

The proof of Proposition 13.11 consists of essentially two steps: first, we prove compactness 
properties in L^(IR^) for the family {V{t))t->q. To that end, we derive uniform bounds with 
respect to r in weighted spaces; this step will be achieved in the current section. Then, we 
prove in the next section, using techniques inherited from dynamical systems theory, that the limit 
of any converging sequence V{Tn) is equal to Fm- As emphasized in the introduction, the proof of 
convergence relies on rather abstract arguments, and thus does not yield any rate of convergence 
in general. However, when the fiux A is linear, the weighted bounds allow us to prove that the 
family U{t) has uniformly bounded moments of order four, and thus ^ holds. As proved in [4], 
the convergence stated in Theorem [1] then takes place with algebraic rate. 

The main result of this section is the following: 

Proposition 3.2. Let C/,™ e L^ n and let U G C([0, oo), ^^(R^)) be the unique solution 

of (fTTj) with initial data C/|i-=o — Uini- Let m > 2{N + 1) be arbitrary, and assume that 

\U,mix)\^{l + \x\^r^^ dx<+00- 

Then there exists a constant Cm > (depending only on m,N, and on the flux A) such that if 
||C^„i||ii < Cm, then 

sup/ \V{t,x)\''{1 + \x\^)'^'^ dx <+oo, 

" rr+1 . (27) 

sup / / \V xV{s, x)\'^ dx ds < +00. 
t>oJt Jr" 

As a consequence, there exists a sequence (t„) of positive numbers such that r„ G [n,n+ 1] for all 
n, and such that the sequence {V{Tn, x))n>o is compact in L^{M.^). 

Moreover, if the flux A is linear, then Cm = +oo for all m > 2{N + 1). 



Before proving the bounds ([271) . we explain how they entail the existence of a converging 
sequence. Thus we admit that (|27|) holds for the time being. First, for any X > 1, r > 0, we have 

\V{t,x)\ dx 

\x\>X 

< (f \V{T,x)f{l + \xfr/^dx] (f il + \x\^)-"^^Ux] 

\J\x\>X ) \.l\x\>X j 

1/2 



sup 

T>0 



|F(t,x)|2(1 + |x|2)"/2 dx^ 
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Since m > TV, we infer that the family {V{T,x)}r>o is equi-integrable. 

Moreover, let /C C be an arbitrary compact set, and let h S be arbitrary, with \h\ < 1. 

Let 

jC:={xe M^, d{x,IC) < 1}. 
The set K. is clearly compact. Then 

\V{T,x + h) - V{t,x)\ dx 



1/2 



< \h\ I [ \VV\{T,x + \h)d\dx 

JK Jo 

< \h\f\VV(T,z)\dz<\h\\JC\'^'(f \W^V{t,x)\ 

Jk. \Jr" 



Now, for all n eN, there exists r„ e [n,n+ 1] such that 

/ |V,y(T„,x)|2dx < / / \V,V{s,x)\'' dxds. 

Consequently, there exists a constant C, depending only on K. and on the bounds on V in 
Lfo(.([0,oo),i/i), such that 

VneN, / |y(T„,a; + /i)-y(T„,a;)|dx< C|/i|. 
Jk 

Hence the sequence {V^(t„, a:)}„>o is equi-continuous in L^(R^). 
Notice also that 

sup||y(T„)|U. <-J—snp\\U{T^)\\L^ < F™^"^ 



n>0 infTN /o n>0 infT" fo 

Thus the sequence {V{Tn,x)}n>o is bounded in L^(M^). 

According to classical results of functional analysis (see for instance [6|), we infer that the 
sequence (F(t„))„>o is compact in L^. 

The rest of the section is devoted to the proof of the bounds l(27|) . We first prove that 
V € L^^{[0, oo), L^{{1 + Then, using the construction of approximate solutions of (flTI) 

performed in the previous section, we derive an energy inequality on the function V. Carefully 
controlling the non- linear terms appearing in this energy inequality, we are led to (|27|) . 

Before addressing the proof, we recall a result which will play a key role in several arguments: 
since Uini G n i^(R^), there exists a positive constant C, depending only on the flux A and 
on ||f7i„i||i, ||[/m.i||oo, such that 

\\U{t)\\l^^9.-, < Ce^^ (28) 



Indeed, performing backwards the parabolic scaling lfTO|) . it turns out that this inequality is equiv- 
alent to the boundedness of u in L°°([0,oo) X M^), where u is the solution of ^ with initial data 
V + Uini- And the bound on u follows from Lemma A.l in the Appendix. 

First step: the family V{t) is locally bounded in L'^{{1 + Ixp)™/^). 

This amounts in fact to proving that U G L^^{[0 , oo) , L"^ {{1 + Hence, multiply (HU 

by U{t,x){1 + and integrate with respect to the variable x. Always with the notation 
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R — e7 , z — Rx + c 
1 d 



R-'-l 
2 ' 



this leads to 



-- / \U{T,x)\'{l+\xn-^'dx 



(29) 



|V,f/p(l + |xn™/2 dx-m {x- W^U)U{1 + dx 
(x ■ V,|C/n(l + |a;n'"/2 _ „^ /" |[/|2|a;|2(l + |a;|2)-i+™/2 rf^. 



B z, 



UiT,x)\ U{t,x) 



RN 



RN 



B z, 



U{r,x) \ _ U{t,x) 
R^ R^ 



■V,U{l + \x\'r/^dy 
• X ?7(1 + |a;|2)-i+™/2 dx. 



Since U{t)/R" is bounded (see l(28|l ). there exists a constant C such that 



B [ Rx + c 



R^-l U{t,x)\ U{t,x) 



R^ 



R^ 



< C 



\Uir,x)\ 
RN ■ 



Moreover, 



c(l + |a;n-l+'"/2 , |x|2(l + |2;n-l+™/2 < (l + |2.|2)m/2 



Hence, using the Cauchy-Schwarz inequahty, we infer that the last two terms in l(29|) are bounded 

by 

]f \V,U\''{l + \x\^r/^ dx + CR^ f \U\''il + \x\T/^ dx. 

4 Jrn Jrn 



On the other hand, 



{x-w,\u\^){i + \x\^r/^dx 

|C/|2 (^N{1 + |a;|2)™/2 + mx^l + \x\^)-^+'''/^^ dx 



< c \uwi + \x\^r/^ dx. 

Gathering all the terms, we obtain 

\U{t, x)\^il + \x\^r/^dx 

N 

< - I \V,U{t,x)\^{1 + \x\T^^ dx + Ce^ 



d_ 



\UiT,x)\''{l+\x\^r/^dx. 



Using Gronwall's Lemma, we deduce that 

U e L,- ([0, oo), L\{1 + \x\^r'^), V^U e LUiO, oo), L2((1 + \x\'r/')). 

Second step: The energy inequality. 

The idea here is the following: assume momentarily that the flux B is Hnear, that is, Bi = 0. 
Let ■0 G i°°([0, oo) X R^) be a solution of (fTTj) such that 'i{}{t, x) > for all t, x. Then, according 
to [17], for any convex function H e C'^(M.), we have 



d_ 
di 



^{t,x)H — 



dx 



H" 



„ ( U (t, x) 

lfj{T, X) 



UiT,x) 



iIj{t,x) 



dx. 
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Taking H : x E M. i-^ x"^ , we infer that 



sup 

T>0 



\U{r,x)\' 



dx 



< +00. 



Hence, if ^{x) behaves like (1 + for |a;| large, the bound in l(27j) is proved. 

Thus, the goal of this step is to build a positive function C/, which behaves like (1 + 
for \x\ large, and which is an approximate solution of the linear part of (fTTj) . with remainder terms 
of order one. Using calculations similar to the ones led in [T^, we then derive an inequality on the 
energy 

U{t,x) ' 



U{t,x) 



U{t,x). 



From now on, we no longer assume that Bi = 0. 

The definition of U is inspired from the construction of an approximate solution in the previous 
paragraph. Precisely, we set 

U{t,x) fQ{z)hm{x) + er^ fi[z) -^/^hmix), T > 0, xeM^, z^e'x + c — - — , 

where the function /i £ T4^i>°°(T^)^ is defined by dH]), and where 

h^{x) := (l + |xn-W2. 

Remember that infTW /o > 0; since 



Vxhm{x) = ~m- 



we deduce that there exists tq > (depending on m), such that 
1 



< -Uz)h^{x) < U{t,x) < 2fo{z)Knix) yyeW\ T> Tq. 



(30) 



We now compute, for r > tq, the rate of growth (or decay) of the energy J \U\'^U ^ . Using equation 
ifTTj) and performing several integrations by parts, we obtain 



d_ 
rf7 



C/(r,x) 



U{t,x) 



U{t, x) dx 



^ U{t,x) 



'U{t,x) 



U (r, x) dx 



UiT,x) 



U{t,x) 
2g(iv+i). 



~drU + A^U + diY^ixU) - e^div:, (("1(2) - c)U^ 



dx 



Bi z, 



U{t,x) 



Va- - ^ ' — ^ dx. 

U{r,x) 



By definition of C/, we have 

-drU + A^U + div^(a;J7) - e^div^ ( (ai(z) - c)U 



div,(x/i„)/o(z) + fo{z)A,h^{x) + 2 ^(z)^'^"^"'^ 



l<i,j<N 



dz 



dxidxj 



[ial,^ - C^)flJ] (z) 



l<i,j<N 

l<i,j<N 



d'^hmjx) 
dxidxj 



—— {XjdxMjrt(X)) + — — 
ox J OXiOX'- 
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where 



z — e X + c- 



Notice that 



divj;(xft,m(a;)) (TV - m)hm{x) + 



(l + |2;|2)i+t' 

and there exists a constant C (depending on m and N) such that for alH, j £ {1, • • • , N}, 



dxidxj 



< C 



(l + |a;|2)i+t' 



\^xhrnix)\ + 



dxidxj 



dxidx^ 



< Chm{x). 



Hence we infer (remember that iV — m < and that inequaUty ((30|) holds) 

-drU + A^U + div^ixU) - e^div^ (("1(2) - c)U^ 

< {N -m)fo{z)h^{x) + Ce-^hm{x)+C- 
N ~m ~ , , „ 1 



1 



(l + |xP) 



2U + ^ 



< 



-U{t,x)+C 



{i + \xvr 



for all T > To, provided tq is chosen sufHciently large. 

On the other hand, since the flux Bi is quadratic near the origin and U je^^ is bounded, we 
have 



,(JV+1)t 



Bi z, 



C/(t,x) 



Nt 



■^J^J^dx 



^ UiT,x) 



■f/(T,x) 



UiT,x) 

dx. 



Gathering all the terms, we obtain 







u 


drj 


IN 


u 



u 



m — N 





u 






2 


/ 




U + 2 f 




u 




u 


JR" 


u 





< c 

+Ce 



dx 



U{;x) 

Ui;x)J (l+NP) 



2U + ^ 



(l-Ar)T 



|C/(.,x)|^ 



da:. 



(31) 
(32) 
(33) 



U{;x) 

Notice that when the flux A is linear, the term l(33|) is zero. 
Third step: control of the term (|32|) . 

Set (/) := U/U ; then according to the first step, <j) G -^i^cd'^o, 00), i^(^m))nL^(,p([To, 00), H^{hm))- 
Moreover, 

V(02/i„) = 2(t>h^V4> + ^^V/i™; 

since |V/im| < m/im, we deduce that cfi^hm G iJo(,([ro, 00), T4^^'^(R^)), and thus, using Sobolev 
embeddings, (j)'^hm e i,i„^([ro, 00), L?' (M^)), where p* := 7V/(A^ - 1) if > 2, and p* = +00 if 
N = 1. Additionally, the following inequality holds: there exists a constant C, depending only on 
TV, such that for all r > tq 



l|0^^m(r)||ip*(HiV) < C||V(0^/l,„(T))||il(RiV) 



< C||0(r)|U.(,„)||V^(r)|U.(,„) + C|l0(T)||i.(, 
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We use the above inequality in order to control the term l(32|) . First, let us write 

|0(r, x)p(l + |xn-(i+^)da; 



where the exponents a, b satisfy 



2a + & = 2 
m ,m m 



which leads to a = 1 — ^, 6 = ^. Notice that a,b £ (0, 1), provided m is large enough (to > 4, 
which is always satisfied if to > 2{N + 1)). 
Then, using Holder's inequality, we infer 

/ \(j){T,x)mi+\x\^)-i^+^) dx 
< ||'^^W'lm||LP(K«) ll'/'M^m|lLl(RiV), 

where the parameter p is given by 



1-^ 

1- A- 



Notice that p is always larger than one. In order to be able to interpolate between and , 
p must also be smaller than p*; if N = 1, p* = oo, and thus we always have p < p* . If iV > 2, 
this condition amounts to to > 2(7V+ 1); we assume that to always satisfies this assumption in the 
sequel. 

Now, let e £ (0, 1) such that 

1 _9 1-e 

pi p* ' 
using once again Holder's inequality, we obtain 

|0(-,x)|2(l + |a:n-(i+^)dx 



+ C| I (/)^/l,n 1 1^1 ||0/im||il(RN-). 

If = 1, then 9 = p^^, and straightforward computations lead to 



a9 + a^—^ = 1- a{l~9)^-. 
2 to to 

Hence, using Young's inequality, we deduce that for all A > 0, there exists a constant Cx such that 

\^iT,x)\^{l + \xf)-i^+'^) dx 



If TV > 2, the calculations are similar and lead to 



1-61 N + 2 , , 2N 
a9 + a = 1 , a{l-9) = . 

2 TO TO 
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Hence l(34| is also valid in this case. 

Using inequality l(30|) and choosing the parameter A small enough leads eventually to 



m < 



m-N 
16 



U{t) 



U{t) 



U{t) 



1 



U{t) 



U{t) 



(35) 



+C / \U{T,x)\dx 



for all r > tq. 

Fourth step: control of the term ([33]). 

Remark 3. We recall that (|33| =0 if the flux A is linear. Hence this step is required only in the 
nonlinear case. 

Using inequality l(28|) . we infer that there exists a constant C such that 



m<c I \u{t,x)\ 



1+^ 



dx. 



U{t,x) 

From now on, we treat the cases iV = 1, TV = 2, and TV > 3 separately, and we set (/) = U /U. 
• If iV = 1, we have, for all r > tq. 



UiT,x) 

V X ~ 



dx 



< C 



\U{t.x)\^+^ . 

C/(t,x) 

(9.0(r)/iV2) (^(^)/,„)i/2 (^^m^^)h^ 



< 



dx = l^'/^h,ndx<l>-cl)^/^dxh„,, 

>x (0^/2/1™) < C||0/l„,||^(' {\\dx<f>\\L^ih^) + Uh^h^)) ■ 



Moreover, 
and thus 

Eventually, we obtain, using once again ([30 

& < C\\U{t)\\li 

Since 



U{t) 



U{t) 



U{t) 



U{r) 



U{r) 



\\U{t)\\l^ <\\U,^,\\li Vt>0 
we infer that if |l(7mi||i is sufficiently small, then 

2 



m < 



m — N 



16 



C/(r) 



U{t) 



Uir) 



U{t). 



(36) 
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If iV = 2, using the Sobolev embedding W^''^{W) C L'^{R^), we obtain, for t > tq 



i[/(.,x)r+^ 



'Ui;x) 



As is the case iV = 1, we have 
Hence we are led to 



< 



dx < C 

< C 

< C 

1/2 



Li(R2) 



c||</>;i™||^/^(||a,0|| 



3l<qiC/(r)|li('(R.) 



U{t) 



U{t) 



U{r) 



U{t) 



U{t) 



Following exactly the same argument as in the case = 1, we deduce that if ||L/i„i||i is 
sufficiently small, then l(36|) holds. 



If TV > 3, we have, for r > tq. 



\U{-,x)\^+- 



yj U{;x) 



< Cl |</)|i+^|V,0|/^ 



'U{-,x) 
1+^ 



dx 



< C 



where the parameter p is such that 



LP(R") 



1 1 1 



1, 



N 2 p 

i.e. p = {2N)/{N ~ 2). Using the Sobolev embedding H^{R^) C LP(M^), we have 



Thus, once again, we obtain 



33D<C||C/(r)|l^/^ 



LP(R") 



< c 



') 



< C(||V,^|U2(,„^) + ||0|U.(,.„)) 



U{t) 



U{t) 



U{t) 



U{t) 



U{r) 



and thus (|36| holds as long as HJ/miiUli is not too large. 



Gathering inequalities (|3T|) . (|35l) and (|36|) . we infer that if ||t/i„i||i is sufficiently small, then for 
all T > To, 



-/ 


{/(r) 




U{t) 



U{t) 



m — N 



/ 


U{r) 






/r" 


U{t) 


U{r) 



< c 



U{t)< 



\U{T,x)\dx] <C||C/„,||ii. (37) 
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Fifth step: Conclusion. 

Let Ci :— (m — N)/4, C2 C\\Uini\\i. Using a Gronwall type argument, we deduce that for 
all T > To, we have 



U{r) 



U{t) 



U{t) 



'Ci(r-s) 



U{s) 



U{s) ds 



Uiro) 



U{ro) 



U{to) + 



< C \U{To.x)\\l + \x\'r" + 



2\m/2 I ^ 



Using l(30l). we infer 



sup 

r>ro 



sup 

r>To 



\u{T,x)Wi + \x\^r'^ <c I \u{T^,x)\\i + \x\^r'^ + c^, 



T+1 



/ 






U{s) 



hmds<C I |C/(to, + |xnW2 + eg 



Hence f7 G L°°([0, 00), ^^((i + |a;|2)W2)), gince /o is bounded away from zero, the bound on V 
follows. 

Concerning the bound on ^ xV , notice that 



V{t,x) 



Uir,x) 
U{t,x) 



hm{x) + e~ 



Mz) 



and thus 



\VxViT,x)\<Ch^{x) 
Consequently, for all t > 



'U{t,x) 



U{t,x) 



U{t,x) 



R" 



\VxViT, x)|2(l + |a;n"/2 dx<C 



U{r) 



U{t) 



which leads to the bound on 'S/xV. (Notice that we even retrieve that VxV € Lf^^ii^, 00), ^^((l + 



Hence Proposition 13.21 is proved. 

Let us now conclude this section by explaining how the bound ([2]) on the moments of order four 
follows from (|27l) . Let Uini S L'^(h^), with m > 2(7V + 2) sufficiently large. Then we have proved 
that U e L°°{[0, 00), _L2(/i~^)). Now, for all t > 0, using a simple Holder inequality, we infer that 



\U{t,x)\\x\'' dx < ||C/(r)||^.(,-i) (^^^ |x|«(l + \x\')-"^/' dx 



1/2 



Hence, if m > iV + 8, we deduce that U e L°°{[0,oo), L^{\x\'^)); going back to the original variables, 
this entails that ^ is satisfied. Thus the convergence result ^ holds if the fiux A is hnear, and 
Proposition [TTT] is proved. 



4 Long-time behaviour 

This section is devoted to the rest of the proof of Theorem [TJ The idea is to use the compactness 
proved in the previous section (see Proposition 13. 2p together with techniques from dynamical 
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systems theory. This type of proof was initiated by S. Osher and J. Ralston in [18], in which 
the authors proved the stabihty of traveUing waves for a quasihnear paraboHc equation. Their 
arguments were then adapted successfully to various kinds of problems in the context of scalar 
conservation laws (see for instance the review in |20j). 

In the present study, our scheme of proof is in fact closely related to the one of M. Escobedo and 
E. Zuazua in [12]; indeed, the idea is to apply the dynamical systems tools to the rescaled parabolic 
system (fTT|) rather than the original conservation law |(T]) . The main difference with [12] lies in the 
presence of highly oscillating coefficients in Ijlip ; thus it is necessary to work simultaneously with 
the homogenized equation (|22| - l|23|) and with the oscillating one. 

Let us now introduce some notation and definitions. First, we denote by St (t > 0) the semi- 
group associated with the homogenized equation, that is equation l(23|) if TV = 1, and equation (|22l) 
if TV > 2. According to Proposition [531 the semi-group Sr is well-defined in L°°(R^) n L'^iK); 
additionally, the contraction property holds, namely 

\\SrFi~SrF2\\Li(M^)<\\Fi~F2\\LimN) Vt > 0, VFi,i^2 e i:°°(M^)nL2(x). 



Hence Sr can be extended on L^(R^). 

We also define the cj-Hmit set associated with a given function Uini <£ L^{R^): recaUing the 
definition of the function V (see l(26|) ). we set 

n[U^m] := {V e L\R^), 3t„ -> c5o, 1/(t„) ^ ^ in L^R^)} , 

where the function U in (|26l) is the unique solution of (fTTj) with initial data Uini- When there is 
no ambiguity, we will simply write Q instead of r2[C/i„i]. 

Notice that V{Tn) converges towards V in if and only if 



lim 

n — ^oo 



U{Tn,x) - fo e'^"a: + c V{x 



dx = 0. 



This equivalence will be used repeatedly throughout the section. 

The organisation of this section is the following: we first introduce a "quasi-Lyapunov function" 
for the semi-group associated with equation ifTTj) . We then prove that Proposition 13.11 holds when 
the initial data Uini has a sufficiently small norm. Eventually, we prove Proposition 13.11 in the 
general case. 

4.1 A quasi-Lyapunov function 

Let us introduce the notion of quasi-Lyapunov function: 

Definition 4.1. Let X be a Banach space, and let H : [0, oo) x X ^ M.. Let {u{t)}t>o be a 
trajectory in X . We say that H is a quasi-Lyapunov function for the trajectory u if the following 
properties hold: 

(i) The family 7i(i, u{t)) (t > 0) is bounded in M; 

(a) There exists a function : [0,oo) — > [0,oo) such that limt_^oo ^{t) = and 

Mt > 0, snp{n{s, u{s)) - n{t, u{t))) < i}{t). 

s>t 

We then have the following result: 

Lemma 4.1. Let X be a Banach space, and let {u{t)}t>o be a trajectory in X. Let Ti : [0, co) xX ^ 
R be a quasi-Lyapunov function for the trajectory u. Then Ti.(t, u{t)) has a finite limit as t oo. 
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Proof. First, since H{t,u{t)) is bounded for t e [0,oo), the quantities 

H := limM n{t,u{t)), H := lim sup it(t)) 

are well-defined and belong to R, with H_< H. 

Let e > arbitrary. There exists > such that 

ip{t) <e yt> te. 

By definition of H_, there exists > such that 

\H{se.u{s,))~H\<e. 

Since 7i is a quasi-Lyapunov function, for all s > s^, we have 

nisM-'^)) < HGs„u(s,)) + v(ss) 

< H + 2e. 

Hence 

H <H + 2e Ve > 0, 
and H = H_. Thus the quantity has a finite Hmit as i — ^ oo. 

□ 

We now apply this notion to the present context: 

Lemma 4.2. Let M e R 6e arbitrary, and let Uim e ^^(R^). For r > and u e ^^(R^), define 
the function Ti. by 

n{T,u):= f \uix)-U^PP[FM]{T,x;en\ dx, 

where the function V^^p was introduced in Definition \2.1[ 

Let U S C(|0, oo), L^(R^)) be the solution of (fTTI) with initial data Ur=o = Umi- 
Then H is a quasi-Lyapunov function for the trajectory {U{t))}t->o in L^{M.^). As a conse- 
quence, the function 

r e [0, oo) I— > / \U{t, x) — fQ{z)FM{x)\ dx, with z — x + c , 

converges as r — > oo. 

Proof. This property is an easy consequence of the first point in Lemma [275} indeed, according to 
Lemma [2751 there exists a constant C, depending only on N and M, such that 

±H{r,U{T)) = ^\\U{r) - U-^^[FM]ir)h^M-) < Ce'^ 

Consequently, for all t' > r > 0, we have 

nir', U{t')) - nir, U{t)) < C{e-^ - e"^') < Ce'^ 
Thus property (ii) of Definition 14.11 is satisfied. Additionally, notice that 

0<7^(t,[/(t)) < ||C/(r)|Ui(R«) + ||/o|Uoo(T«)||FM|lLi(R") 
+ Ce"^ (||VFm||li(r") + ll-P'M|li2(Riv)^ 

+ Ce"^^ [\\Fm\\w^,1{RN) + ||FM|li3(RiV) + ||VFAf||^2(RiV)) 

< l|C^jm||Li(B«) + C- 
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Whence 'H{t,U{t)) is bounded for r € [0,oo). Consequently H is a quasi-Lyapunov function 
for the trajectory U{t). According to Lemma [4.1^ H{t,U{t)) admits a finite limit as t — > oo. 
Furthermore, we have 

||;7(t) - fo(z)FMhr - \\U{t) - C/^PP[i^M](T) + e-^Ui + e^'^C/zH^i , 

where Ui and U2 are defined by (fr7l) - (|20l) and (pT]) respectively. Hence for all r > 0, there holds 

H(t, C/(t)) ~ Ce-^ < ||(7(r) - fo{z)FMhi < n{T, U{t)) + Ce-^ 

where the constant C depends only on W^''^ bounds on Fm- Thus the function 

r i-> ||[/(r) - /o(^;)i^M||ii(Rjv) 

converges as t ^ 00, and 

lim \\U{t)~ f^{z)FM\\L^ = lim H(t,C/(t)). 

□ 

Definition 4.2. Lei Uini G i^(R^) 6e arbitrary, and let M :— Uini. Let U he the solution of 
ifTTj) wii/i initial data U\t^Q = Uini- 
We define the number i{Uini) by 



lim 

r — >oo 



\U{t,x) — fQ{z)FM{x)\ dx, with z — e^x + c- 



=2t 



Notice that Proposition [3?T] is equivalent to 

^{U^ni) = Vt/„„; G ^^(R^). 

Classically, we now derive a continuity property for the function t. 
Lemma 4.3. The function 

U e Li(R^) ^ e{U) e M 

is Lipschitz continuous. 

Proof Let ?7£],i7^] e L^iR^), and let M^*) = for i = 1,2. We denote by U^'^ G 

C([0, 00), Li(M^)) the solution of ^ with initial data U^^^. Then for all r > 0, the contraction 
principle ensures that 



C/(l)(^)_{/(2)(^) 

Hence, for all r > 0, we have 

/ U^'Hr,^)- fo{z)FM^,j{x) 

< 



< 



Li(R") 



da; 



U^^\t,x)^ fa{z)FMi2-, {x) 



Li(K") 



dx 



Li(K") 

According to Lemma A. 2 in the Appendix, 

||i^A./(i) - Fj,^(2) llii(Riv) 



< 
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Eventually, we obtain, for all r > 0, 



f/(i)(r,x)-/o(z)FM(i)(a;) 



dx - / (r, x) - fo{z)FM(V (x) dx 

< (l + ll/o|lioo(T«)) 

and thus, passing to the limit. 
Hence ^ is a Lipschitz continuous function. 



< 



im ini 



L1(R«) 



Li(K™) 



□ 



4.2 Analysis of the tu-limit set 
Proposition 4.1. Lei Ui„i G ^^(M^), and sei 

M := 



Assume that the uj-limit set associated with Umi is non-empty. Then the following properties 
hold: 

(i) For all V C,VL, 



V = M; 



(ii) SM C n for all t > 0; 
(iii) For all F G il, we have 



Proof. Throughout the proof, we denote by U the unique solution of equation (fTTI) with initial 

Property (i) is quite straightforward: indeed, conservation of mass for the equation (fTTI) imphes 
that 

/ U{t) =M Vt > 0. 
If F G fi, then there exists a sequence (t„)„>o such that 

lim T„ = oo and lim / \U{Tm x) — fo (z„) V"(a;) dx — 0, 



where Zn — e'^'^x + c- 



According to a result of G. Allaire (see [2]), 



hm / /o (z„) V(x) dx = (/o) / V; 



gathering the three equalities, we obtain property (i). 

We now address the proof of property (ii) , which rehes on the second point in Lemma 12. 5j let 
y G f7 be arbitrary, and for all e > 0, let G L'^{K) n L°°(R^) such that 



\Ve - "^^||li(R«) < e. 
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Let (r„)„>o be a sequence of positive numbers such that r„ ^ oo and 

\U{Tr,,x)- fo{Zn)V{x)\ dx^O, 



where z„ = e'^"-x + c 



by 



Let p e C^(R ) be a mollyfing kernel; for S > 0, set ps := S " p{-/S), and define the function 



Then Lemma 12.51 ensures that U^''^ satisfies equation (fTTj) with an error term, the latter being 
bounded for all T > in L°°([0, T], L\R'^)) by 

UJT,e{S) + CT,s,Se~'^" 

where uJT,e '■ [0, oo) [0,oo) is such that limo+ ujt,s — 0, and where the constant CT,e,s depends 
only on e, 6, N and T. 

Using the contraction principle for scalar conservation laws, we infer that for all T > 0, and 
for all T £ [0,T], 



dx 



\U{Tn+T,x)-U^f{T,x)\ dx 



< lutA^) + CT,e.se'^" + / \U{t,,,x) - Jo{z^)V{x)\ dx 



u'Sr^^{x) - h{zn)V{x) 



dx. 



Now, according to Definition 12.1 



U^;'{t,x) = {SrV,)*^ps{x)h[e^-+^x + c 



g2(r„+r) _ I 



,-(r„+r). 



U\ T, X, e^" X + c 



„2(r„+r) _ I 



Hence for all t e [0, T], we have 



2 I r, x, e'^"^'^x + c 



g2(r„+r) _ I 



U^f{T,x)-SrVh\e^-+^x + c 



„2(r„+T) _ I 



dx 



< ll/olloo sup llS'^y - *x P5||Li(RiV) 

re[0,T] 

+e^'^" (I|C^i||l = ([0,T]xTJ',L1(R^)) + l|t^2||L°=([0,T]xT«,Li(R^))) 



< ll/oll 



sup ||5rV;-(5.v;)*. 

Te[o,T] 



+e~'^" (II'^i||l = ([0,T]xTJ',L1(R^)) + ll^^2||L~([0,T]xTf ,Li(R^))) 

Gathering the two inequalities, we deduce that for all r S [0, T], for all n, 5, e, 

p2(r„+r) _ 



U{t^ + t, x) - S'r1^(a;)/o e^"+^x + c 



< UJT,e{5)+CT,e.Se + \U (jn, x) - fo{Zn)V {x)\ dx + Ce. 
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In the right-hand side of the above inequaUty, we first choose e sufficiently small, then 6 so that 
ujT,e{5) is sufficiently small, and eventually n large enough so that the two remaining terms are 
small as well; hence 

lim inf (cjT.eW + CT,e,5e"^" + / \U{Tn,x)- h{zn)V{x)\ dx + Ce\=Q. 

n-^co £>0,(5>0 V Ion I 



Thus we have proved that for all T > 0, 



lim sup 



C/(t„ + t, x) - (5.y)/o e^"+^x + c 



j2(r„+r) _ Y 



dx = 0. 



The above convergence entails immediately that SrV G SI for all t G [0,T]. Since T > was 
arbitrary, property (ii) is proved. 

There remains to prove property (iii), which is a variant of the LaSalle invariance principle; let 
V G Q.he arbitrary, and let r„ be a sequence of positive numbers such that lim„^oo t„ = +00 and 



lim 

n — *oo 



\U{Tn,x) ~ fo{z„)V{x)\ dx = 0, 



where z„ = e'^"x + c^-^ — -. According to a result of G. Allaire (see [2]), we have, since (/o) = 1 
and /o e C(T^), 



V-F, 



M\ 



lim 

n — *oo 



fo{Zn)\V{x)-FM{x)\dx 

= lim / foizn) 



Vix)~^^j^ + ^^j^-FM{x) 



h{Zn) fo{Zn) 



dx 



= lim / fo{zn) 



U{t„,x) 



fo{z„ 



Fm{x) 



dx 



Consequently, 



Corollary 4.1. Let U^m e Li(R^), and set 



□ 



M := 



Assume that the u-limit set n[Uini] is non-empty. Then £(Uini) — 0, and thus the result of 
Proposition lg.il holds. 



0; 



Proof. Let V G fl he arbitrary. Then 

lim II^.^-FmIIli 



this property is stated in Proposition 12.11 in the case when V G L°°(M^) n L^{K), but can be in 
fact easily generalized to an arbitrary function V G hy using the contractivity of the semi-group 
Sr-. indeed, let e > 0, and let K & L°°{R^) n L^iK) such that JV^ = JV = M, and 

IlK — ^||li(r«) < £■ 

Then for all r > 0, 

WSrV^FMh < \\SrV ~ SrZWi + \\SrV, ~ FmWi < \\V ^ + \\SrVe - Fm\\i. 
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Hence, using Proposition 12. 11 we infer that 



\imsnp\\SrV - Fm\\i <e Ve > 0, 

r — *oo 

and thus \\StV — Fm\\i vanishes as t cx). 

On the other hand, property (ii) in Proposition 14.11 ensures that StV £ for all r > 0, and 
thus, using (iii), 

\\SrV - FmWi ^ iiU„u) Vr>0. 



Consequently, £{Uini) — 0. Going back to the definition oi £{Uini), we deduce that 

dx = 0. 

□ 



lim 

T—^OO 



C/(t, x)-fo[ ex + c^-—^ ] Fm{x) 



Thus the proof of Proposition 13. II is complete provided we are able to show that the set ^\Uini] 
is non-empty for a sufficiently large class of functions Uini € L^(M^). In the case when ||?7mi||i is 
small, this result follows from Proposition 13.21 and from a contraction principle. The proof in the 
general case is more involved, and in fact, an analysis similar to the one performed in Section [3] 
has to be conducted once more. 

4.3 Proof of Proposition [3TT] when ||?7i„i||i is small 

We now complete the proof of Theorem [T] when is small. Let Uini S L^(M^). Assume that 

Uini satisfies the following assumptions 

3m>2(A^+l), ||C/„,||ii(Kiv) < C„,, (38) 
and U,m e LH(1 + \x\^r/^) n L°-(E^), (39) 

where the constant Cm was introduced in Proposition 13.21 Then according to Proposition 13. 2^ the 
w-limit set r2[[/mi] is non-empty, and consequently Proposition 13.11 is true (see Corollarv l4.ip . 

Let us now prove that Proposition [3lT] holds when Uini merely satisfies (|38|) : this fact is a direct 
consequence of the density of L'^{{1 + |a:;|2)™/2) n L°°(R^) in Li(R^), together with the continuity 
of I. Indeed, for all £ > 0, let C/f„ G ^^((l + |a;|2)™/2) n L°°(R^) such that 

\\Uini - C^i>ii||Li(R") < l|f^fmll-Li(K") < CJm- 

Then £{Uf^^) — 0. Since £ is Lipschitz continuous (see Lemma rO)) . there exists a constant C such 
that 

£mm) = m^m) - £{UU\ < C\\U,m ^ UUIl^M-) < Cs. 



Since the above inequality holds for all e > 0, we deduce that £{Uini) — 0. RecaUing the definition 
oi £, we infer that Proposition 13.11 holds for all initial data Uini G i^(M^) satisfying fM)) . 

4.4 Proof of Proposition 13.11 in the general case 

The case when ||J7i„i||i is large follows from the following Lemma: 

Lemma 4.4. There exists a constant Cq, depending only on N and on the flux A, such that for 
all U^nt G L^K^), 

£{U,m) <Co^ mm) = 0. 
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Before proving the above Lemma, let us explain why the result of Proposition 13.11 follows. If 
t{Uini) < Co for all Uini S i^(R^), then the above Lemma states that I is identically zero on 
L^(R^), and thus Proposition 13.11 is true. Thus we assume by contradiction that there exists 
Uini G i^(K.^) such that £{Uini) > Cq. Consider the function 

(f>:t& [0,1]<^ £{tU^n^). 

We have proved in the previous paragraph that (j) vanishes in a neighbourhood of zero. Moreover, 
(/) is a continuous function according to Lemma rOl Now, it is obvious that 0(0) = 0, and 0(1) = 
(■{Uini) > Cq. Hence there exists to G (0, 1) such that 

Hto) = ^ 

But according to Lemma [T4l 0(to) — i{toUini) = 0, which is impossible. Thus ({Uini) — for all 

U^n^ G L\^^ ) . 

There remains to prove Lemma 14.41 According to Corollary 14.11 and using by now standard 
arguments, we only have to prove that there exists a set C L^(M^), which is dense in L^{W^), 
and such that 

3C > 0, yUrm e l\r^) n A, £{u,m) <c^ n[u^m] ^ 0. (40) 

In the following, we will take A — L^{{1 + \x\'^)"^^^), for some m > sufficiently large. 

The scheme of proof of the implication l(40|) is very similar to the one of Proposition [321 indeed, 
we have to prove that if £{Uini) is small enough, then there exists a sequence (r„) of positive 
numbers, with lim„^(x) t„ = +oo, such that (F(r„, •))„^oo is a compact sequence in L^(R^). 
Notice that this is obviously equivalent to the compactness of the sequence F(t„, •) — Fa/, whose 

norm is of the order of i{Uini) as n ^ oo. Thus our strategy is the following: rather than 
using directly the equation on U, we consider the equation on the function U — J7^PP[FAf]. We 
prove that for an appropriate function U, an inequality of the type l(3T|) holds, with U replaced by 
U — U^^^[Fm]- Then, all the occurrences of |1C/(t)|1i in the proof of Proposition 13 . 21 are replaced by 
II ([/— J7*pp[Fa/])(t)|| 1, which converges towards £{Uini) as r ^ oo. Thus the same arguments which 
led us to compactness in the case when ||C/mi||i is small show that compactness holds, provided 
£{Uini) is small enough. 

Let us now retrace the main lines of the proof: first, consider a function Uini G L^(M^) such 
that Uini G L^((l + |a;p)™/^) for some sufficiently large m (to be chosen later). Set M = J^^, Uini 
and 

WiT,x) = U{t,x) - C/^PP[FM](T,a;;e^). 

In the rest of the proof, for the sake of brevity, we will write [/'^pp (r, a;) as a short-hand for 
J7^pp[Fm](t, x; e"^). Then the following properties hold 

W G L,- ([0, ex.), L'{il + |xn™/2)) n LL([0, ex.), ((1 + IxH'-/^)), 
3C > 0, Vr > 0, ||W^(t, OIIl-cr") < Ce^\ 
lim \\W{T)\\Li(RN)=£{U,m). 

Moreover, using Lemma [231 we deduce that W satisfies 

drW = diy^{xW) + A^w-Rdiy,j:{{ai{z)-c)W) 



-i?^+Miv, 



B, U, V (z, ^^"^ 



with R = , z ^ Rx + ^ , and we recall that the remainder JJ^*^"^ satisfies 



|C/™'"(t)|Uoo(k«) + ||C/™"(t)||^.(,„.,.) < Ce-- (41) 
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for some 7 > 0. 

Then, using the bounds on U, C/^pp together with the regularity assumptions on B, it can be 
easily proved that 



Bi z, 



U{t,x) 



Bi z, 



2a2(z)/o(z) -rrj + 5(r, x), 



J12N 



and the function b is such that there exists C > such that 



V(T,a;) e M+ X M", \b{T,x)\<C 
We define a function W by 



W{t, x) 



R 



-2N-1 



\W{t, x)\ 



W{t,x) ^Wo{x,z) + e~^Wi{x,z), 
with Wo{x, z) ~ fo{z)h„i{x) and 

-A,Wi + div,(aiW'i) = 2Ay,Wo - div,((ai - c)Wo) - ljv=i2div^(a2/oi^A/Wo). 
Notice that by definition of /o and c, the compatibility condition is always satisfied, and 

Wi{x, z) = /i(z) • Vyhra(x) + 1n=iWi{z)Fm{x)Ki{x) , 

with 



Let To > such that 



-A^wi +<l\Y^{aiWi) = -2div2(a2/o)- 



W{t, x) > ^fo{z)h„,{x) Vr > To, Vy G M^. 



For further purposes, we also choose tq such that 

|1M^(t,-)||i <2^([/„„) Vr>ro. 

(Notice that if £{Uini) = there is nothing to prove). 

Using calculations similar to the ones performed in the proof of Proposition I3.2[ we infer that 
for T > To, 



-/ 


w 


dr Jrjv 


w 


^/ ( 


W{ 







W{t) 



m — N 



dx 



W 
W 



\WiT,x)J (l + |a;|2)i+t 



+Ce 



{1-N)t 



\W{t,x)\^' 



V 



W{t,x) 



W{t) + 2 



dx 



W 



W{t) 



1 


W[t, x) 


/r« 


W{t,x) 



|C/'''^"(r,a;)| dx. 



Using the same arguments as in the third step of the proof of Proposition 13. 2| we deduce that if 

m>2{N + \), 



W{t,x) 



R« \W{t,x)J (l + |a;|2) 



dx 



< 



m- N 
16 

1 

+ 2 



W{t) 



W{t) 
2 



W{t) 



W{t) 



-C£{Urmf 
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Similarly, the calculations of the fourth step in the proof of Proposition 13.21 yield 

W{t,x) 



e(i-^)- / \W{t,x)\ 



/ 


W{t) 




W{t) 



""Wir^x) 



dx 



W{t) 



W{t) 



W{t) 



W{t) 



W{t) 



W{t) 



W{t) 



Eventually, using the Cauchy-Schwarz inequality together with the bound (|4T|, we infer that 

|i7'^'=™(T,x)| dx 

Wir) ' 



/ 


Wir, x) 




W{t, x) 



< C\\U'^'"'{t)\\ 


L2(e-'l = l^ 




W{t) ' 







W{t) 
Wir) 



1/2 



< C 



m- N 
16 



/ 


W{t) 




W{t) 



W{t) 

\ 1/2 

W{t)\ 



W(t). 



W{t) 



W{t) 



1/2 



Gathering all the terms, we deduce that there exists a constant Cm, depending only on N and 
m, such that if £(Uini) < Cm, then for all r > tq. 



-/ 


W{t) 


dr Ju« 


W{t) 



W{t) 



m — N 
16 



Wi7 



W{7 



W{t) 



R« 



W{t) 



W{t) < C. 



Compactness of a subsequence VF(r„) follows. Hence the w-limit set is non-empty, and thus 

mrn) = 0. 



Appendix A 

Lemma A.l. Assume that the flux A satisfies ([5]). Let v G W^'°° {T'^) be a periodic stationary 
solution of (51), and let u G L^^([0, oo), L°°(R^)) n C([0, c»), Li^^(K^)) he the unique solution of 
^ with initial data mm G v{y) + n Then u G L°°([0,oo) x R^). 

Proof. This result was proved in [7] in the case iV = 1. When N > 2, the proof goes along the 
same lines; the only difference lies in the use of the Sobolev embeddings, which depend on the 
dimension. Hence we merely recall here the main steps of the proof, with an emphasis on the case 
N>2. 

In the rest of the proof, we set f{t, y) — u{t, y) — v{y). Then / solves the equation 

dtf + divyB{yJ)-Ayf^O, (42) 

and according to ^ the flux B is such that for all / G M, 

|div,B(2;,/)| <C(|/| + |/r), 
\dfB{yJ)\<C{\f\ + \fn, 
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where the exponent n is such that n < {N + 2) /N. Moreover, 

II/(^)IIl1(R") < ll"mi - w||l1(R") > 0. 

For q>l arbitrary, multiply (|42|) by |/|'', and integrate over R^. Using a few integrations by parts 
(see [7j), we are led to 



d 
dt 



I/I 



9+1 



v,|/|' 



< Co 



I/I 



9+1 



I/I 



q+n 



(43) 



We then use Sobolev embeddings in order to control the L'^^^ and L*+" norms in the right-hand 
side. We distinguish between the cases N — 2 and > 3, since the space is critical in dimension 
two. 

• If TV 2, then H^R"^) C ^^(R^) for all p G [2,oo). Interpolating between and LP 

for some p sufficiently large, we have 



II/IIl.+"(r^) < ll/ll?ll/llp-' with 



1 



1 - 



q + n 1 



< 11/11? i/r 



< c,\\m vi/i 



_2p_ 

51+1 

9 + 1 



Notice that 



and 



q + n 



(i-e) 



q + n — I 



'^ + 1" ' (9 + l)(l-i)' 
q + n— 1 



9+1 



< 1 > 1 



since n < {N + 2)/N. Thus, we choose p > 1 such that 

q + n — 1 



(<? + !) (l-^) 



< 1. 



Young's inequality then implies that for all A > 0, there exists a constant C\,q and an exponent qi 
such that 



i/r+" <A vi/|V +cxjf\\r 



(44) 



The other term in the right-hand side of l(43l) can be bounded in a similar fashion: we have, for all 
A > 0, 



1/1^+1 < A V|/|V +CA,,||/||f, 

R« 2 



(45) 



for some exponent q2 which can be explicitely computed. Choosing an appropriate parameter A, 
we infer that there exist qi,q2 > such that 



d 
di 



r" 



I/I 



9+1 



V,|/|V <C,(||/||f -i-ii/iif). 



Using l(45|) one more time leads to 



dt 



-rJ ifr'+cj i/r^<Q(ii/iif + ii/iif)<c. 



R" 
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Using a Gronwall-type argument, we infer that / G L°°([0, oo), L's (R^)) for all q> 1. 
• When TV > 3, we use the Sobolev embedding i/i(R^) C Lp' (M^), where 

2N 

P' 



iV-2 



Interpolating L^^'^ between and L 2 ^ we obtain 



<ii/ii^(,+i)ii/iir' 


< 


I/I = 


2£) 

g+1 


11/11} 










P* 


2£) 




< c 


VI/ 


■j+i 
2 


g+1 1 



l/ll}- 



where the parameter 9 £ (0, 1) is given by 

1 



26 



1 - I 



It can be checked that 



q + n p*{q + I) 

N + 2 e{q + n) 
N ^ q+1 



n < 



< 1. 



Hence (|44|) holds when > 3. Inequality l(45l) is proved with similar arguments. As in the two- 
dimensional case, we deduce that / £ L°°{[0, oo), L^(R^)) for all q. Using Theorem 8.1 in Chapter 
III of [T6] (see [7J for details), we infer eventually that / £ L°°{[0,oo) 



□ 



Appendix B 

Lemma A.2. Let M > M' be arbitrary. Then 

FM{y)>FM'{y) VyeR^. 

As a consequence, 

\\Fm-Fm'\\i = M-M'. 



Proof. The arguments are exactly the ones which lead to the uniqueness of stationary solutions of 
l(22|) . l(23|l . and they can be found in [l]. We recall the main steps below for the reader's convenience. 
Let F Fm - Fm' ■ Then F £ (1 C2(R^), and F > 0. Hence the set 

e := {x £ R'^,F{x) > 0} 

is non-empty. The idea is to prove that F+ — FIq satisfies a Hnear elliptic equation; since > 0, 
F+ cannot vanish anywhere, and thus F+{x) > for all x £ R^. 

Let us now derive an equation on F+ . Substracting the equations on Fm and Fm' , we have 



E 



d^F 



dxidxj 

l<i,j<N ' ■> 



div,(&F) =0, 



where 



b{x)^a{FMix)+FM'ix))-x, x £R^; 
notice that a = if TV > 2. Since F £ H'^{R^), we have 

div:,(6i^)le -divx(6F+) 
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almost everywhere. Thus, we obtain 

l<i.j<N ' ^ 

Integrating the above equation on leads to 



dxjdxi 



^ 0. 



Let us now perform the change of variables (|25| . which changes the matrix rj into identity: setting 
F{y) = F{Py), and 9 := {F > 0}, we infer 



[AyF = c[ 



0. 



dxjdxi 



® l<i,j<N * 

Moreover, F e H'^ n W'^''^{R'^), and thus Lemma 7 in [I] applies. We deduce that 

AyiFl^) = leA.F, 

and thus 

l<i,j<N •' l<iJ<N ■' 

Eventually, F+ solves the elliptic equation 

l<i,j<N ' ^ 

with h G LJ^^(R^). Using either a unique continuation principle or Harnack's inequality (see |13) . 
Theorem 8.20), we infer that if F+ vanishes at some point x in M/^ , then F+ is identically zero on 
R^, which is absurd. Hence F+{x) > for all x S M^, and thus \ 9 = 0, which means that 
F{x) > for all x e R^. 

□ 
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